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Abstract —The problem of three-dimensional (3-D) wind profile 
prediction is addressed based a trinion wind model, which 
inherently reckons the coupling of the three perpendicular 
components of a wind field. The augmented trinion statistics are 
developed and employed to enhance the prediction performance 
due to its full exploitation of the second-order statistics. The 
proposed trinion domain processing can be regarded as a more 
compact version of the existing quaternion-valued approach, with 
a lower computational complexity. Simulations based on recorded 
wind data are provided to demonstrate the effectiveness of the 
proposed methods. 

Index Terms —wind profile prediction, trlnion-valued repre¬ 
sentation, least mean squares, adaptive filtering, widely linear 
processing 

I. Introduction 

Hypercomlex numbers are high-dimensional extensions of 
real numbers ID, 0, and they have been introduced to solve 
multivariate signal processing problems, such as colour image 
processing 0-161, vector-sensor array signal processing 171- 
ca, human gesture spotting ca, ca, wind profile predic¬ 
tion ITSl - llTl . l22l . and wireless communications ll2a . In 
the last application, anemometer readings are modeled with 
hypercomplex numbers and the wind profile is predicted by 
adaptive filtering algorithms. In particular, pure quaternions 
have been widely used to model three-dimensional (3D) wind 
speed. When external atmospheric parameters are available, a 
full quaternion-valued model can be considered US). 

The quaternion-valued model leads to improved perfor¬ 
mance over real-valued models, since it accounts for the 
coupling of the wind measurements and can be extended to 
exploiting the augmented quaternion statistics ifTSl . However, 
pure quaternions do not form a mathematical ring IfTSl . since 
the product of two pure quaternions is not a pure quaternion 
in general. As a result, the related adaptive algorithms for 3-D 
wind profile prediction initialised in the pure quaternion will 
have to work in the full quaternion domain. The prediction 
results have to be truncated from full quaternions to pure 
quaternions, implying redundant computations in the update 
process. Eor example, it takes 16 real-valued multiplications 
and 12 real-valued additions to implement a full-quaternion¬ 



Fig. 1. Trinion-valued three-dimensional wind speed model. 

valued multiplication. By comparison, it only takes 9 real¬ 
valued multiplications and 6 real-valued additions to calculate 
the multiplication of any two 3-D numbers. 

In this work, we aim to develop a more compact 3-D wind 
speed model based on a 3-D mathematical ring called trinions, 
termed by an anonymous author who provided all possible 
definitions of 3-D numbers 03 . We first define the gradient 
operation in the trinion domain and then derive the least mean 
squares (LMS) adaptive algorithm, which is further extended 
to augmented statistics. 

The rest of this paper is organised as follows. A brief 
introduction to trinions is provided in Section II. The trinion- 
valued LMS (TLMS) algorithm is developed in Section III, 
and the augmented TLMS (ATLMS) algorithm in Section IV. 
Simulation results are presented in Section V and conclusions 
are drawn in Section VI. 

II. Trinions 

As shown in Eig. 1, the 3-D wind speed is a tri-variate 
signal composed of three perpendicular components, and can 
be modeled with 3-D hypercomplex numbers (see Eig. 1). 
There are various definitions of a 3-D number v composed 
of one real part {va) and two imaginary parts {vb and Vc), 

V = Va + IVb + JVc, (1) 

according to definitions about the relationships among the 
three base elements 1, i, and j. 




To form an Abelian group of these three, the following rules 
apply 1201: 

= = = -^■ (2) 


Trinions subject to the above rules form a commutative ring, 
namely, for two trinions v and w, we have vw = wv. 

The modulus of v is given by l20l 

kl = (3) 

We define the following conjugate of v, 

V* =Va- JVb - IVc, (4) 


SO that jup is equal to the real part of vv*, i.e. |up = 

To develop an LMS-like algorithm in the trinion domain, 
the trinion-valued gradient needs to be defined. In the complex 
domain, a variable z and its conjugate can be viewed as two 
independent variables, based on which the complex-valued 
gradient can be defined l2TI . Similarly, in the quaternion 
domain, a variable q and its three involutions can be viewed 
as four independent variables, based on which the quaternion¬ 
valued gradient can be derived Il22ll . l24l . However, to our best 
knowledge, the trinion involution is not available in general. In 
this paper, we define the following gradients of a function fix) 
with respect to the trinion-valued variable x and its conjugate. 


V./= 

V.. / = i /-f zV. J-f jV. J), 


(5) 


respectively, where x = Xa + iXb + jx^. 

Unlike the quaternion-valued gradient, the trinion-valued 
gradient remains the same, no matter which side of the sub¬ 
gradients the imaginary units i and j are on, since trinions are 
commutative. We can then calculate the derivatives of some 
simple functions, for instance. 


dx _dx* 
dx dx* 

dx dx* 1 — i + j 

dx* dx 3 

III. Trinion-valued LMS Filtering 


( 6 ) 


Similar to the LMS algorithm in the complex-valued do¬ 
main, the trinion-valued error e{n) is given by 

e(n) = d{n) — w'^{n)x{n), (7) 

where d{n) is the reference signal, w{n) is the weight vector, 
x[n) = [x{n — P)\x{n — P — 1); • ■ • ; x{n — L — P + 1)]^ is 
the filter input, L is the filter length, and P is the prediction 
step. The cost function is expressed as 

J(n) = |e(n)p. (8) 


Using the steepest descent method, the following gradient need 
to be calculated 

V^„.J(n) = ■^[V^„„J(r^) +iV^u^J{n) +jVw^J{n)]. (9) 


For details of the calculation, please refer to Appendix at the 
end. Using results there, the update of the weight vector can 
be obtained as 


win + 1) = w{n) + iie{n)x*{n), (10) 

where /i is the step size. The algorithm formulated above is 
termed as the trinion-valued LMS (TLMS) algorithm. 


IV. Augmented Trinion Statistics 


A zero-mean trinion-valued vector v is composed of three 
zero-mean real-valued random vectors Va, Vb, and Vc, and 
their complete second-order statistics can be found in the 
following six real-valued covariance matrices: 

E{VcvJ}, = E{VavJ}, ( 11 ) 

= E{V]jV^ }, C'llolla ~ }• 

And they can be more efficiently represented by three trinion- 
valued covariance matrices: 

= E{vv^}, 

C,,. = E{vv^^}, (12) 

= E{vv^^}, 

where superscript ^ denotes the Hermitian transpose and the 
two additional mappings of v are defined as 

V" = Vb - IVa - JVc, 

= Vc-lVb- JVa- 


Neither of these two mappings is an involution, and they are 
defined as shorthand notations only. Then we can obtain 


^ VaVa 




= 

,z jC VV^ ) 

^VcVc 


'U'lyd ) , 

^ Va^b 


jC ) , 

CvbVc 

= 

1 J^VV^')-! 



'iC VV^- 


To take the complete second-order statistics into consideration, 
we need to make use of the augmented input vector composed 
of the original input vector a;(n) and its two mappings 
x'^{n),x^in), namely. 




xin) 
x\n) , 
x^{n) 


(15) 


and the predicted estimate yin) will be given by 
i/(n) = 

= wjin)xin) + w'^in)x^in) + wjin)x\n), 


(16) 







Fig. 2. Predicted result from the TLMS algorithm. 


where = [wi-,W 2 ;w 3 ]. Analogous to the derivation of 

the TLMS algorithm, the update equation of the augmented 
TLMS (ATLMS) algorithm is given by 

+ 1) = (17) 

where p is the step size. 

The computational complexity for each update of the weight 
vector of the trinion-based and quaternion-based filtering algo¬ 
rithms are summarised in Table 1. It can be seen that the trinion 
model can effectively reduce the computational complexity 
compared to the quaternion model. 


TABLE I 

COMPUTIONAL COMPLEXITY PER UPDATE OF THE WEIGHT VECTOR 


Algorithm 

Real multiplications 

Real additions 

QLMS 1171 

16L-I-4 

16L 

Augmented QLMS 

64L-I-4 

64L 

TLMS 

9L-I-3 

9L 

Augmented TLMS 

27L-I-3 

27L 


V. Simulations 

In this section, both proposed algorithms (TLMS and 
ATLMS) are applied to anemometer readings provided by 
Google’s RE<C Initiative ||25]| . The wind speed measured on 
May 25, 2011 is used for demonstration. The step size is set 
to be 6 X 10“^. The filter length is 8, and the prediction step 
is 1. All algorithms are initialised with zeros. The predicted 
results provided by TLMS and ATLMS algorithms are shown 
in Figs. 2 and 3, respectively. From the results, we can see 
that both algorithms can track the wind data effectively. 

The learning curves averaged over 200 trials of the proposed 
algorithms are shown in Fig. 4, compared with the quaternion- 
based QLMS and AQLMS algorithms. It can be observed that 
both augmented algorithms (AQLMS and ATLMS) have a 



Fig. 3. Predicted result from the ATLMS algorithm. 



Fig. 4. Averaged learning curves. 


faster convergence rate than the original ones (QLMS and 
TLMS), due to their full exploitation of the second-order 
statistics. Meanwhile, the proposed TLMS algorithm has a 
similar performance with the QLMS algorithm, while the 
ATLMS algorithm is comparable with the AQLMS algorithm. 
However, as shown in Table I, the proposed trinion-based 
algorithms have a much lower computational complexity. 

VI. Conclusion 

A compact model for three-dimensional wind profile pre¬ 
diction based on trinion algebra has been proposed, with two 
LMS-type adaptive filtering algorithms derived using partial 
and full trinion-valued second-order statistics, respectively. 
Numerical simulations using recorded wind data have shown 
that the two algorithms have a similar performance to their 
































































quaternion-valued counterparts, but with a much lower com¬ 
putational complexity. 


Appendix 

The cost function J (n) as a function of real-valued variables 
is given by 


J = {da- wjxa + w'^Xa + Xh)'^ 

+ {db - wjxb - wjxa + wjxc)"^ (18) 

-f {dc - wjxc - wjxb - w'^Xa)'^, 

where the time index ‘n’ has been dropped for the sake of 
compact notation. Then the three component-wise gradients 
can be computed as 


J = 

2[{XaxJ + XbxJ + X, 

T\ 

+ 

{xbxl 

+ XcxJ 

- XaxJ)Wb 

+ 

[XcX^ 

- XaxJ 

- Xbx]:)wc 

- 

[daXa 

+ dbXb + deXe)] 

1 , 

J = 

2[{XaxJ + Xbx]: - Xc 

T\ 

+ 

f T 

{XcX^ 

+ XaX^ 

+ Xbx'^)Wb 

+ 

[XcIP'h 

T 

- XaX^ 

+ XbxJ)Wc; 

+ 

{daXc 

— dbXd — 

- deXb)] 

1 

J = 

2[{XaX^ - XbxJ - Xc 

T\ 

,x^ )w, 

+ 

{XbX^ 

T 

+ Xax\ 

'i)wb 

+ 

( T 

+ XcxJ 

+ XbxJ)Wc 

+ 

{da^b 

"h dbXc, ~ 

- deXa)] 



(19) 


( 20 ) 


( 21 ) 


Finally, we can obtain the expression for the gradient of J(n) 
by substituting (fT^- dSlT i into @, 

J(n) = ^e(n)a;*(n), (22) 

which yields the update equation in ( fTOb . as the term | can 
be absorbed into the step size. 
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